Novel Phase Between Band and Mott Insulators in Two Dimensions 
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We investigate the ground state phase diagram of the half-filled repulsive Hubbard model in two 
dimensions in the presence of a staggered potential A, the so-called ionic Hubbard model, using 
cluster dynamical mean field theory. We find that for large Coulomb repulsion, U ^> A, the system 
is a Mott insulator (MI). For weak to intermediate values of A, on decreasing U , the Mott gap closes 
at a critical value Z7 c i (A) beyond which a correlated insulating phase with possible bond order (BO) 
is found. Further, this phase undergoes a first-order transition to a band insulator (BI) at C/ C 2(A) 
with a finite charge gap at the transition. For large A, there is a direct first-order transition from 
a MI to a BI with a single metallic point at the phase boundary. 

PACS numbers: 71.10.-w, 71.10.Fd, 71.10.Hf, 71.27. +a, 71.30.+h, 71.45.Lr 
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The Mott metal-insulator transition is a paradigmatic 
problem in condensed matter physics that has been real- 
ized in a broad set of experiments, ranging from the high 
temperature superconductors to recently observed phe- 
nomena involving cold atoms trapped in optical lattices. 
Mott insulators (MI) are characterized by a density of 
one electron per unit-cell and strong Coulomb interac- 
tions that often leads to symmetry breaking via antifer- 
romagnetic (AF) order. In contrast, charge localization 
in band insulators (BI) occurs when the number of elec- 
trons per unit cell is even, so that the bands are either full 
or empty. These two types of insulators can be realized 
within a single system, the Ionic Hubbard model (IHM) 
at half-filling, by tuning the strength of Coulomb inter- 
action. The BI-MI quantum phase transition in the IHM 
has been well studied in the limit of either oneij2iM or 
infinite dimensions^, but significant uncertainties remain 
in the phase diagram. Originally suggested to describe 
charge-transfer organic salts^, there has been renewed 
interest in the model due to potential applications to fer- 
roelectric perovskites 7 -. Confinement of Fermionic atoms 
in an optical lattice to realize the ionic Hubbard model 
is also conceivable in the future, because recent experi- 
mental progress has allowed for Fermionic atoms to be 
cooled well below the degeneracy temperature^. 

Our objective in this letter is to obtain the ground 
state phase diagram of the IHM at half-filling in two 
dimensions using cluster dynamical mean field theory 
(CDMFT) 9 . The IHM is realized on a bipartite lattice 
by the assignment of an alternating chemical potential 
(±A) for each sublattice on top of the regular Hubbard 
model. The Hamiltonian can be written as 
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Here, t and U denote the hopping amplitude between 
nearest-neighbor sites and the onsite Coulomb repulsion, 
respectively. +A(— A) denotes the local potential energy 
for the A(B) sublattice. We are interested in the case 



where the combined density of A and B sublattices is 1 
and so we set the chemical potential, /i, to {7/2. In the 
atomic limit (t = 0) and for {7/2 < A, the ground state 
has two electrons on the B sublattice and none on the A 
sublattice resulting in charge density wave (CDW) order 
with a band gap of A — {7/2. In the opposite limit where 
{7/2 > A, each site is occupied by one electron and a MI 
is formed with a gap U. Therefore, in the atomic limit at 
exactly U c = 2A the system is gaplessi 2 -^. This tran- 
sition from a BI to a MI is likely to persist when t ^ 0, 
but non-trivial charge fluctuations in the intermediate to 
strong coupling regime ({7 ~ 2A) render static mean-field 
and perturbative treatments of the transition invalid. 

Bosonization calculations for the IHM in one 
dimension 4 , valid in the weak coupling regime (A -C U <C 
t), find a transition from a BI to a bond-ordered (BO) 
insulator at U = U c , where the charge gap vanishes. The 
BO phase is characterized by finite charge and spin gaps. 
Increasing U further, the spin gap vanishes at U = U s 
in a Kosterlitz-Thouless-type transition leading to a MI. 
Numerical studies for long chains using the density ma- 
trix renormalization group (DMRG) 2 '- in the strong cou- 
pling regime find that the charge gap remains finite at the 
transition point U c to the BO phase, although the optical 
gap does vanish. DMRG is unable to identify a second 
transition point towards a MI because extrapolation to 
large system sizes in the critical region is a challenge to 
numerical methods. 

Our key findings for the IHM in two dimensions are 
the following. A nonzero staggered chemical potential 
induces charge density wave (CDW) order over the entire 
range of {7, with the CDW order parameter approaching 
zero asymptotically for large U. In the large U regime, 
the system behaves as a MI with AF order. For weak 
to intermediate A, decreasing the value of {7, the system 
undergoes a transition to a correlated insulating phase 
suggesting BO with a closing of the charge gap at the 
transition point {7 c i(A). The transition also manifests 
itself via an abrupt increase in the slope of the staggered 
density, double occupancy, and staggered magnetization 
as a function of U. Further, the BO phase undergoes a 
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first-order transition to a BI at a lower critical coupling 
[7 C 2(A) where the charge gap remains finite. AF order is 
found only in the MI and BO phases. For A ~ 4.5t, the 
system undergoes a direct first-order transition from the 
BI to the MI, with a finite charge gap everywhere. 

CDMFT is a non-perturbative technique where the 
full many-body problem is mapped onto local degrees 
of freedom treated exactly within a finite cluster that 
is embedded in a self-consistent bath. It is a natural 
generalization of single-site DMFT to incorporate spa- 
tial correlations. The method has passed rigorous tests 
for the ID Hubbard model where it compares well to ex- 
act solutions^. CDMFT has also been used to elucidate 
various aspects of the phase diagram of the cuprates, 
including the pseudogapii and superconducting phases 
within the 2D Hubbard models. 

Using CDMFT, the IHM on the infinite lattice in two 
dimensions reduces to the cluster-bath Hamiltonian be- 
low, that is subject to a self-consistency condition: 

)• (2) 

■ma mfj,<T 

/i, v = 1, • • •, N c denote indices labelling the cluster sites 
and m — 1,- • -,Nb represent those in the bath. The 
self-consistent calculation proceeds by an initial guess for 
the cluster-bath hybridization V mfla and the bath site 
energies e ma to obtain the cluster Green's function G^ v . 
Applying the Dyson's equation, S c = G^ 1 — G~ x , where 
Gq c denotes the non-interacting Green's function for the 
Hamiltonian in Eq. [2l the cluster self-energy is obtained. 
S c is then used in the self-consistency condition below to 
determine a new Gq c : 



G£{z)-Y. {z) 



Nc_ 
4tt 2 



dK- 



1 



■H - t(K) - E c (z) 



(3) 

Here, K denotes a momentum vector in the reduced Bril- 
louin zone of the cluster superlattice and z = iu) n is the 
Fermionic Mastubara frequency. From Go c in Eq. [3] a 
new set of V m ^a and e mcr 's is generated closing an it- 
erative loop. In this work, the Lanczos method is used 
to solve the cluster-bath Hamiltonian. The cluster size 
is fixed to N c = 4 and the bath size to Nb — 8. The 
Lanczos method can access both the strong and weak 
coupling regimes with equal ease and it is well suited to 
compute dynamical quantities directly in real frequency. 
Rotational symmetries of the cluster on a square lattice, 
together with particle-hole symmetry at half-filling, re- 
duce the number of bath parameters significantly. These 
bath parameters can depend on spin, allowing for sym- 
metry breaking solutions. For details of the method, we 
refer to earlier work2ii£. 

We analyze our results by starting with Fig. [TJ which 
shows the density difference between sublattices A and 
B as a function of U, for A—2t. Clearly, the system 




FIG. 1: Staggered charge density as a function of U in the 
BI, MI, and BO regions for A = 2t. Inset shows an enlarged 
view of the abrupt change of slope at the MI-BO boundary. 
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FIG. 2: (a) Staggered magnetization, and (b) double occu- 
pancy as a function of U, for A = 2t. 



presents long-range charge density wave (CDW) order 
over the entire range of U. The CDW order parameter 
asymptotically approaches zero for increasing U but re- 
mains finite everywhere. Starting from the MI phase at 
large U, the staggered charge density increases smoothly 
as U is reduced, but at U=U c i=5At there is an abrupt 
increase in the slope. This anomaly is clearly identified 
in the inset and it is believed to be a signature of the 
transition to a Bond-Ordered (BO) phase. This result is 
similar to observations in earlier DMRG calculations for 
the one dimensional IHM 2-3 . The BO phase extends up 
to [7={/ C 2=4.6t, where there is a first-order transition to 
a BI signalled by a jump in the staggered density. Start- 
ing from U = 0, the BI shows a continous decrease in 
the staggered density and persists up to U=6.55t. Con- 
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sequently, we have the coexistence of a BI and BO for 
4.6* < U < 5 At and a BI and MI for 5 At < U < 6.55*. 
We computed the total energy in the coexistence region 
(not shown) and find that the BI solution for the CDMFT 
equations is always higher in energy than the BO or the 
MI phase for the same U. 




FIG. 3: Local density of states as a function of U for A = It 
shown with a broadening factor of 0.02t. 

To further characterize the nature of the transition let 
us consider the staggered magnetization, M=n^ — 
and the double occupancy, D=(n^ni), as a function of U 
for A=2t shown in Fig. [2] With reducing U both these 
quantities show an abrupt increase in slope at the same 
transition point U=U c i=5At (see inset) where the BO 
insulator emerges from the MI. Both M and D present a 
jump at the first-order transition from the BO phase to 
a BI. The staggered magnetization abruptly goes to zero 
upon entering the the BI phase signifying an absence of 
AF order. The double occupancy increases monotoni- 
cally with reducing U in the BI state due to the increas- 
ing tendency towards CDW order. 

In Fig. [3l results for the local density of states (LDOS) 
as a function of U for A = 2t are presented. Note that 
the LDOS shows particle-hole symmetry consistent with 
half-filling because it is averaged over the A and B sub- 
lattices. For U equal to the bandwidth of 8*, the spectra 
show lower and upper Hubbard bands consistent with a 
MI. In addition, sharp coherent peaks are observed (ab- 
sent in single-site DMFT). They correspond to the mo- 
tion of a carrier in an AF background and define the Mott 
gap. As U is lowered, the Mott gap reduces its value until 
it closes at the transition point U — 5 At, which mani- 
fested itself earlier as an abrupt change in slope for the 
staggered density, double occupancy, and staggered mag- 
netization. With further reduction in U, the trend for the 
Mott gap is reversed, with an increase in the gap in the 
BO phase. Further, at U = 4.55 there is an abrupt jump 



in the gap as compared to U = 4.60 signalling a tran- 
sition to a BI, which also shows an increasing gap with 
decreasing U . For A > 4.5t, the BI-MI transition has no 
intervening phase and it is accompanied by the closing of 
the charge gap. This transition scenario is reminiscent of 
the phase diagram of the extended Hubbard model in one 
dimension, where an intermediate phase with bond order 
is stabilized between the CDW and SDW phases only for 
intermediate strength of Coulomb repulsion U when the 
nearest neighbor Coulomb repulsion V is varied^. 




FIG. 4: Momentum dependent density-of-states along the 
(0, 0) to (-7T, 71") direction for (a) the BI, (b) the phase boundary 
of the MI-BO transition, and (c) the MI for A = It. 

To characterize the distinct phases observed, we have 
computed the momentum-dependent density-of-states 
using the following relation between the the lattice 
Green's function and the cluster self-energy: 



G(k, z) 



N c ^ 



ik-(jLi— u) 



i 



z + fjt- f (k) - S(z) 



(4) 

In Fig. IU we show the evolution of the density of states, 
A(k, u>), in the first quadrant of the Brillouin zone, from 
(0,0) to (7r, 7r). At large U (Fig. Efc)) sharp coherent 
peaks that define the Mott gap are observed. In addi- 
tion, the lower and upper Hubbard bands are present, as 
well as high energy features separated from them by A. 
The minimum of the gap is found at k=(7r/2, tt/2) where 
the low energy features disperse towards zero frequency, 
while the maximum of the gap is found at k=(7r,7r) or 
(0, 0). For small U (Fig.QJa)), where we have a BI, a dis- 
persion characteristic of AF order is observed but with 
the notable absence of Hubbard bands. By contrast, for 
intermediate U, on the phase boundary of the MI-BO 
transition at U = 5.4*, we found a metallic density of 
states with distinct higher energy bands at (ir, 0) and 
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(0, 7r) that emerge from spontaneous dimerization and 
persist into the BO phase. Therefore, the BO phase, al- 
though similar to the BI in the behavior of the charge gap 
as a function of U, shows a distinct signature of strong 
correlations. The metallic peak itself splits to open a gap 
right below the phase boundary. 
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FIG. 5: Momentum dependent density-of-states along the 
(0, tt) to (tt, 0) direction for A = 2t. 




FIG. 6: Proposed CDMFT phase diagram for the 2D IHM. 

In Fig. A(k,uj) is shown, from (0, tt) to (tt, 0). The 
uncorrelated BI at U — 0M (Fig.[5ja)) shows no disper- 
sion of the low lying features and the gap becomes in- 
dependent of momentum. For the MI phase at U — lOt 
(Fig. Oc)) we see an inward dispersion in the occupied 
band at k=(-7r/2, tt/2) where the gap is a minimum, con- 
sistent with commensurate AF fluctuations. In contrast 
to this result, Fig.[5](b) for U = 5.2t shows a maximum of 
the charge gap at k=(7r/2, tt/2) while the minimum lies 



at k=(0, 7r) or {w, 0), consistent with the spontaneously 
dimcrizcd BO phase with incommensurate order. 

To summarize, we have studied the quantum phase 
transition between a BI and a MI in the two dimensional 
IHM using CDMFT. CDW order persists throughout the 
entire range of U. For A < A c ~ 4.5i, the large U Mott 
insulator state undergoes a transition to a correlated in- 
sulating phase BO with reducing U, accompanied by the 
closing of the charge gap at the transition point t/ c i(A). 
The BO phase undergoes a first-order transition to a 
BI at t7 c2 (A) upon further reduction of U. As shown 
in the phase diagram depicted in Fig. [6j the width of 
the BO phase shrinks as both small and large A are ap- 
proached. For A > A c , a direct first-order MI-BI tran- 
sition with no intervening BO phase is observed. AF 
order is found only in the MI and BO phases. Our re- 
sults differ significantly from earlier ones that used the 
single-site DMFT technique^ where a paramagnetic so- 
lution must be enforced even though the ground state is 
AF. Single-site DMFT misses spatial correlations crucial 
to the physics of this model. Moreover, the single-site 
DMFT finds a finite metallic region between the BI and 
MI whereas we find a single metallic point at the MI- 
BO (A < A c ) or MI-BI (A > A c ) transitions. Longer 
range interactions could give rise to a metallic phase be- 
tween the insulators^. Our results are fairly similar to 
those obtained for this model using DMRG 2 ' 3 in one di- 
mension. An important difference is that at the metallic 
point we observe a closing of the charge gap whereas 
DMRG finds only a closing of the optical gap, with the 
charge gap remaining finite over the entire range of U. 
The two-transition scenario moving from the BI to the 
MI was observed earlier in the weak coupling regime us- 
ing Bosonization^ for this model in one dimension. We 
have shown that a similar transition persists in two di- 
mensions, even in the strong coupling regime when A is 
comparable to half the bandwidth. Modulation of the 
local potential with an alternating value could be real- 
ized in an two dimensional optical lattice loaded with 
ultracold Fermionic atoms. Noise correlations in time-of- 
flight images can be used to distinguish between the AF, 
BI and BO phases^. Energy gaps can be obtained, in 
principle, using rf spectroscopy^. We expect our results, 
predicting a novel correlated insulating phase with bond 
order, to motivate such experiments to verify our finding. 
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